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ESTIMATION OF STICE-FIXED NEUTRAL POINTS OF AIRPLANES 

By Maurice D. White 


SUMMARY 


A method is given for calculating the stick-fixed 
neutral point of an airplane with propeller windmilling, 
flaps neutral, and landing gear retracted. This method 
differs from those formerly used principally in the 
procedure for estimating the effect of the windmilling 
propeller. Comparison of the neutral points predicted 
by this method with neutral points obtained in flight 
tests indicates good agreement at low lift coefficients. 
The methods presented. In conjunction with the results 
given in NACA CB No.. lij.H01 ’’Effect of Power on the Stick- 
Fixed Neutral Points of Several Single- Engine Monoplanes 
as Determined in Flight," should be useful in estimating 
the 3tick-flxed neutral points of new designs for all 
flight ponditions. 


INTRODUCTION 


Since the publication of reference 1, in which a 
method was presented for predicting the static longitu- 
dinal stability of airplanes, additional flight data on 
the longitudinal stability of airplanes have become 
available. Attempts to correlate these longitudinal- 
stability data with the longitudinal- stability data com- 
puted on the basis of reference 1 Indicated that the 
methods of reference 1 were inadequate when applied to 
unconventional designs. A more rational method has 
therefore been developed for computing the longitudinal 
stability of airplanes In terms of the stick-fixed neutral 
point, which yields results In good agreement with flight 
results. This method differs from the method of refer- 
ence 1 chiefly in the procedure for estimating the effects 
on longitudinal stability of the windmilling propeller 
and of the fuselage and nacelles, although the differences 
in change in longitudinal stability due to the fuselage 
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and nacelles as computed by the two methods is generally 
small for conventional designs. The methods given may 
also he adapted readily to the estimation of the effect 
on neutral point of changes in configuration of existing 
designs. 


SYMBOLS 


a 0 two-dimensional lift -curve slope, per degree 

w width of fuselage or nacelle at selected longitu- 

dinal station, feet 

c wing chord, feet 

o' mean aerodynamic chord, feet 

D propeller diameter, fe6t 

/l 

q dynamic pressure, pounds per square foot ( — -pV c y 

r factor used in correcting lift -curve slope for 

effect of end plates 

p air density, slugs per cubic foot 

V velocity, feet per second 

S gross area, including section through fuselage, 

square feet 

M pitching moment, foot-pounds 

x distance along longitudinal axis, feet 

xj_ distance from wing leading edge to middle of 
* fuselage section, feet 

Xj distance from wing leading edge to front o.f 

fuselage section directly -ahead of wing leading 
edge, feet 

x 0 neutral-point location, percent c 

Ij. longitudinal distance from center of gravity to 

quarter-chord point of mean chord of horizontal 

tail, feat (J^CLASSIFIEO 



NACA OB NO. L5C01 


13 


a 

£ 


Cm 

°N 

A 


longitudinal distance from center of gravity to 
propeller plane 

angle of inclination of air flow relative to 
longitudinal axis, radians 


angle of attack relative to free stream, radians 
angle of downwash, radians 
lift coefficient 


/ Li ft\ 

Ws / 


pitching-moment coefficient 
normal-force coefficient 


C 


Vise ; 

Normal force 


qS 


■) 


aspect ratio 

thrust disk-leading coefficient 



T thrust, pounds 

N number of propellers., nacelles, or fuselages 

Subscripts: 

w wing 

f fuselage 

nao nacelle 

t horizontal tail 

p propeller 

p-j_ propeller normal force 

P 2 downwash due to propeller 

L.B. wing leading edge 

mid .midchord point of local chord 

T.E. wing trailing edge 



k 
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COMPUTATION OP NEUTRAL POINT 


The neutral point of an airplane is defined as the 
center -of -gravity location at which the slope of the curve 
of airplane pitching-moment coefficient against lift coef- 
ficient dCm/dC^ is zero* This airplane pitching-moment 
slope is the resultant of the .pitching-moment slopes 
contributed by the various : parts of the airplane, in 
order to predict the neutral point of an airplane, the 
estimated values of. dCm/dCL due to the various parts of 
the airplane may be combined at each of several center- 
of-gravity locations and the neutral point may then be 
established as the center-of-gravity location at which 
the resultant value of dC^/dCj, is zero. A procedure 
equivalent to this procedure for determining the neutral 
point is Illustrated In figure 1, where the values 
of dCju/dCk for the various parts of the airplane 
are plotted against center-of-gravity location. The 
center-of-gravity location at which the total of the 
positive values of dCm/dCi, is equal to the negative of 
the total of the negative values of dPm/dCL - that is, 
at which dC m /dCT is equal to zero - is the neutral 
point. 

Detailed procedures for calculating the values 
of due to the principal parts of the airplane - 

wing, fuselage, nacelles, horizontal tail, and propellers - 
are given herein. The methods apply best at low lift 
coefficients where the airplane drag coefficient may be 
considered not to vary with lift coefficient, the 
parameters involved In the calculations are most nearly 
linear, and the effects of air-flow separation are at a 
minimum. The airplane was assumed to operate with pro- 
peller windmilling, flaps neutral, and landing gear 
retracted. Negative changes in pitching -moment slope are 
stabilizing and positive changes are destabilizing. 

n’he values of £C m /dCL for individual parts 
computed by the following methods are based on the values 
of dCj/da for the wing alone. In order to determine 
the actual value of dCm/dCi, based on the total lift of 
the airplane the values of dCfc/dCj, computed by the 
present methods should be multiplied by the ratio 
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This correction does not affect the location of the 
neutral point as determined by the present methods because 
at that center-of-gravity location the resultant value 
of d Cm/d C l of zero would be unaffected by the correction. 

The methods described in the present report are 
effective for predicting stick-fixed neutral points at 
low lift coefficients with propeller windmilling, flaps 
neutral, and - landing gear retracted. Inasmuch as the 
stick-fixed neutral point generally remains fixed or 
moves back with increasing lift coefficient in this con- 
ditioh of flight (reference 2), the results obtained from 
the present methods will be conservative for the higher 
lift coefficients. 

The effects of power on the neutral-point location 
cannot at the present time be predicted by methods com- 
parable to those given herein. By the use of data given 
in reference 2, however, it should be possible to make 
reasonable preliminary estimates of the shift in stick- 
fixed neutral point due to power. 

A knowledge of the elevator hinge-moment character- 
istics is essential for determining the shift in neutral, 
point due to freeing the elevator* 

Pi tching-moment slope due to wing .- At a given 
center-of -gravity location the value of the pitching- 
moment slope due to the wing with flaps neutral 
(dQm/6CL) w is numerically equal to the distance between 
the center of gravity and tl.e wing aerodynamic center 
expressed as a fraction of the mean aerodynamic chord. 

For the present purposes It has been found - satisfactory 
to consider that the wing aerodynamic center Is located 
at a percentage of the mean aerodynamic chord equal to 
the average of the percentages of the chords at which the 
aerodynamic centers of the root and tip airfoil sections 
are located. For center-of -gravity locations behind the 
wing aerodynamic center, (dCm/dCL^ is positive. 
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Pitching-moment slope due to fuselage and to engine 
n acelles .- The value of due to the fuselage or 
to the nacelles may he expressed by 



fl) 


The value of 

methods 
value of 


l/dM\ 

q(4, 


of reference 
1/dMN 


m)f. 


nac 


is calculated directly by the 

nac 

3. According to this method the 
may be computed from 



(2) 


where the integral is taken along the entire length of 
the fuselage or nacelle. 

It Is convenient in evaluating equation (2) to 
divide the fuselage or nacelle into finite sections* 

The factor Is then calculated for each section 

2 da 

with average values of w and dp/da used for each 

section and, finally, the values of ^w^^Ax for all 
the sections are totaled. . 

The factor dp/da represents the variation with 
angle of attack of air-flow angle relative to the X-axis 
of the fuselage or nacelle; in the calculation of dp/da, 
it is assumed that fuselage - or nacelle - interference 
effects may b'e ignored. Between the wing leading edge 
and trailing edge,. where the flow follows the wing surface 
dp./da is considered to be zero. Behind the wing trailing 
edge the value of dp/da is assumed to vary linearly from 

zero at the wing trailing edge to ^1 - at t* 16 tail. 

The discrepancy between the assumed linear variation and 
the actual variation of dp/da behind the wing will 
normally be of little importance because, for conventional 


NACA C3 No. L5C01 


7 


fuselage arrangements, this portion of the computed 
fuselage moment is only about 10 percent of the entire 
fuselage moment. 


Ahead of the wing, values of dp/da are greater than 
1.0 because of the upwash induced by the wing; the variation 
of dp/da with distance from the wing leading edge is 
given in figure 2, which is taken directly from reference 3 % 
Figure 2(b)* is used for all sections except the section 
directly ahead of the wing leading edge where' values of 
dp/da rise sharply as the leading edge is approached. 

For this region an average integrated value of dp/da. 


defined as- 


3s 

da 


used, 
of x 


- 1 .r 1 ^ 
da 


dx and given by figure 2(a), is 

In determining cLp/da in figure 2(e.) the value 
used in calculating the abscissa is the distance 


from the wing leading edge to the front of the .fuselage 
section; in determining dp/da in figure 2(b) the value 
of x-i -used in calculating the abscissa is the distance 

from the wing leading edge to the middle of the fuselage 
section considered. 


The curves in figure 2 correspond to a wing lift- 
curve slope of lj-5 per radian. In order to correct for 
other values of the lift-curve slope the values given by 
the curves are increased or reduced in direct ratio to 
the lift- curve' slope. 

An additional factor accounting for the effect of 
the fuselage or nacelle on the wing becomos important 
when the width of the fuselage or nacelle is not uniform 
along the wing chord. According to reference J this 
factor is computed by ‘ 

| HE = N lJ( w L.B. + . 2w mid - 5 W T. E . )° 2 (5) 


When the plan form of the body shows local protu- 
berances In width directly ahead of the wing leading edge 
the application of equations (2) and (J) Is believed to 
result in exaggerated values of the destabilizing effect. 
For such cases, reduction of the increment due to the 
protuberance i obtained by calculating the factors with 
and without the protuberance) by one-half appears to 
result In more nearly correct values t This correction, 
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which is "based on very scant data, is subject to modifi- 
cation as more data are obtained. 

Curves for estimating the wing lift -curve 
slope dCT/da as a function of aspect ratio are given in 
figure 5. These curves were obtained from the expression 

ao L _ a 0 * 

da , . ao x S7 .5 

I + r 

tiA 


end at • low aspect ratios were corrected on the basis of 
data given in reference I 4 , The correction for the effect 
of vertical tails mounted as end plates ffig. I}.) was a ^ so 
taken from reference I 4 . . 

The values of &Q m /'bCi J due to the fuselages and 
nacelles are always positive. 


Pitc hing-moment slope due to hor?.zontal tail .- The 
value of the pitching-moment slope duo to the horizontal 
tail ( 6 C m / 6 CL)^ is calculated by the expression 



q t 

“q da fc 


St*t(l - 
7 .£Cl 

S * c da 


* 1 ) 

day t 


(h) 


The various terms in equation (if) are calculated as 
follows j 


±t 

q"* 

dCu^. 
<3 at 


A value of 0,9 is used for the condition of 
wi ndm i 1 1 i ng pr ope Her . 


Curves for the evaluation of 
figures 3 a nd if.* 


«3C W . 
Hb 

flat 


are riven in 


2.1 Values of e at the tail may be estimated from 

fla the charts of reference 5 , A shorter empirical 

method for determining a c/da at the tail, 
based on these charts, is given in reference 1. 
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For conventional tail locations the value of dCm/d'OL 
due to the horizontal tail is negative. 

Pitching -moment slope due to normal force on 
windmllling propeller and to downwash increment resulting 
from this normal forc'e .- The values of d.G^/dCjj due to 
the propeller are considered to arise from two sources - 
the normal force on the inclined propellerjand the 
resulting increment of dovmwash at the tail. The value 
of dOm/dCx, due to the normal force on the propeller is 
given by 



The terms in equation (5) are calculated as follows i 


( 


dG iA 

da/ p 


Reference 6 gives methods for estimating the 
/ dCw\ . 

value of ( J that should he used when the 
geometric and operating characteristics of the 


propeller are known. (In reference 6 [ - 


is designated Cy 




.) 


\ ^“/p 


Unpublished flight data 


obtained at the Langley Laboratory on a number 
of airplanes indicates that with propeller 
windmilling and at the airspeeds for which the 
present methods are applicable the propeller 
thrust coefficient T c is on the average 
approximately -0.03 . The data also show that 
when the propeller pitch setting is set for 
high rotational speeds the propeller blades 
are, at these airspeeds, still against the low- 
pitch stops, which correspond on the average 
to a propeller-blade angle at O.75 propeller 
radius of about 20°. It has been found suffi- 
ciently accurate for the airplanes considered 
in the present analysis, to use directly the 


values of 


( 

\ TT-J -for a propeller blade angle 
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of 20° given in figure 2 of reference 6, This 
figure presents data for a Hamilton Standard 
J155-6 propeller at a value of T c equal to 
zero and a side-force factor (defined in 
reference 6) of 80.7. For ‘convenience the 
values given in this figure for a propeller 
blade angle of 20° are tabulated as follows*. 


Propeller 

\ da / p 

Two-blade 

0.095 

Three -blade 

.155 

Four -blade 

.170 

Six -blade (single rotating) 

.240 

Six -blade (dual rotating) 

.295 


It should be emphasized, however, that when 
the propeller characteristics are known the 
proper values of side-force factor and blade 
angle should be used in reference 6 to deter- 
mine the value of 




Values of are obtained from figure 2(b), 

\aa/p 

A value of that corresponds to the plane . 

of the propeller is used. 


For propeller locations ahead of the center of 
gravity the value of dC m /6CL due the propeller normal 
force is positive. 


The added downwash over the horizontal tail resulting 
from the normal force on the propeller contributes a 
''urther change in pitching -moment slope, that is assumed 
0 be given by 



(0) 
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Equation (6) is derived by calculating the downwash in 
the slipstream from momentum considerations and reducing 
the value by an empirical’ coefficient to take into account 
the effects of the wing, fuselage, and other factors. In 
all cases in which a 1 value of 6C m /dC;g exists because of 
this factor - that is, when the horizontal tail may be 
considered as being within the downwash field due to the 
propeller, as on conventional tractor -propeller arrange- 
ments with tail at rear of fuselage - the value will be 
positive. 


The increased rate of change of downwash over the tail- 
accounted for by equation (6) should also be considered in 
estimations of the absolute angles of attack of the hori- 
zontal tail or of absolute elevator trim angles. In order 
to include this effect, the rate of .change of horizontal 
tail angle of attack with v/ing angle. 'of attack 

da t _ /. de\ 
do“ " - 3a) t , 

should be reduced by the quantity 



The change in absolute downwash angle due to this effect 
will of course bo zero when the resultant angle of attack 
of the propeller axis, including the effect of wing 
upwash, is zero. 


The effect on the bodies within the influence of 
the propeller wake of the added downwash due to the pro- 
peller was ignored. 


COMPARISON WITH PLIGHT DATA 


The methods described in the preceding section and 
illustrated in figure 1 have been applied to the estima- 
tion of the stick-fixed neutral points of the airplanes 
shown In figure .5; dimensions of these airplanes are 
given in table I. The estimated values of dCm/dCx, due 
to each component and the estimated resultant neutral 
point are given for each airplane in -table II. The stick- 
fixed neutral point for.fa'ch airplane as determined from 

w— w - UNCLASSIFIED 
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flight tests for a value of C^ = 0.3 using the method 
described in reference 2 is also shown in table II. 

Plight values and estimated values of neutral-point 
location are compared in figure 6 . The results generally 
agree within ± 1.5 percent c, which represents the accu- 
racy with which the flight neutral points shown are 
believed to have been evaluated. The lack of agreement 
between estimated and flight neutral points for airplane 9 
may be due to the fact, shown by wind-tunnel tests (refer- 
ence 7 ), that at normal flight attitudes the horizontal 
tail of this airplane is within an unusually low-energy 
wing wake. 


CONCLUDING REMARKS 


The present report gives a method of estimating the 
stick-fixed neutral point of an airplane with propeller 
windmilling, flaps neutral, and landing gear retracted. 

The method may be applied best at low lift coefficients. 

A study of available flight data indicates that at higher 
lift coefficients the neutral point generally remains the 
same or moves back; thus, the neutral point predicted by 
the present method will be conservative throughout the 
flight range for this flight condition. 

No reliable methods for predicting the effects of 
power ane at the present time available. For preliminary 
estimates of the shift in stick-fixed neutral point due 
to power, however, the data given in NACA CB No. L4K01 
"Effect of Power on the Stick-Fixed Neutral Points of 
Several Single-Engine Monoplanes as Determined in Flight" 
may be applied. 

The shift in neutral point due to freeing the elevator 
may be computed if the elevator hinge-moment character- 
istics are known. 


Langley Memorial Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va. 
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TABIE I.- DIMENSIONS OF AIRPLANE TESTED 
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Air- 

plane 

Test- 

point 

symbol 

Wing 

area, 

% 

(sq ft) 

Wing 

span 

(ft) 

Aspect 
ratio 
of wing, 
,A w 

Mean 
aero- 
dynamic 
chord 
of wing, 

(S 

(ft) 

Hori- 

zontal 

tall 

area, 

St 

(sq ft) 

Hori- 

zontal 

tall 

span 

(ft) 

Aspect 
ratio 
of tall 

A) 

Longi- 
tudinal 
distance 
from wing 
a.c. to - 
1/4 mean 
chord of 
horizontal 
tall 
(ft) 

Pro- 

peller 

dia- 

meter, 

D 

(ft) 

Longi- 
tudinal 
distance 
from wing 
a.c. to 
propeller 
plane 
(ft) 

Fuse- 

lage 

length 

(ft) 

Maximum 

fuselage 

width 

(ft) 

. 

Nacelle 

length 

(ft) 

Max- 

imum 

nacelle 

width 

(ft) 


O 


34-0 

7.5 

4-75 

28.0 

9.3 

3-1 

11.97 

6.00 

6.04 


3-5 

— 

— 


+ 


57-3 

5-9 

6.80 

48.0 

12.8 

3.4 

16.47 

10.00 

6.71 


4-1 

— 

— 


X 

236 

37-3 

5-9 

6.80 

48.0 

12.8 

3-4 

16.48 

11.08 

8.93 

31.7 

3-2 

— 

— 


a 

213 

34.0 

5-4 - 

6.72 

40.0 

13.0 

4.2 

15-52 

10.38 

9.65 

30.6 

2.8 

— 

— 


o 

300 

40.8 

5-5 

7-29 

58.3 

16.0 

4-4 

19-83 

12.17 

8.71 

35-3 

4-5 

— 

— 


A 

236 

37.0 

5-8 

6.64 

4 i -4 

13.2 

4.2 

15-34 

10.76 

9.10 

32.2 

2.8 

— 

— 


V 

. 2h& 

38.3 

5-9 

6.88 

45-2 

13-3 

3-9 

15-44 

11.58 

10.84 

32.9 

2.9 

— 

— 

8 

V 

328 

52.0 

8.3 

7.02 

54.0 

21.8 

8.8 

21.46 

11.50 

7-72 ' 

19-5 

5.1 

33.8 

3-5 

9 

N 

258 

39-3 

6.0 

6.94 

61.2 

14-8 

3.6 

15-79 

9.00 

7.56 

26.9 

4-9 

1 

— 

10 

IS. 

33 k 

42.8 

5-5 

8.11 

77-8 

18.5 


19-43 

13.08 

8.50 

33-8 

5-1 


— 

11 

A 

422 

49-7 

5-9 

9-11 

107.4 

19.O 

3-4 

17.90 

' 12.00 

11.04 

36.7 

4-7 


— 

12 

> 

490 

54-2 

6.0 

9-76 

110.8 

20.8 

3-9 

22.69 

13.00 

U.56 

41.0 

5.0 


— 

15 

< 

. 959 

89.5 

8.4 

11.74 

203.0 

26.7 

3-5 

32.40 

11.50 

11.93 

56.O 

6.3 

15-4 

4-3 

14 

X ) 

46 5 

61.3 

8.1 

8.36 

100.0 

21.2 

4-5 

24.02 

11.25 

5.81 

47.9 

4-3 

19.2 

4-0 

15 


538 

61.3 


9.68 

122.3 

20.2 

3-3 

25.58 

12.08 

7-48 

48.5 

3-5 

16.9 

4-1 

16 

y> 

. 540 

70.0 

II 

8.13 

116.1 

23.1 

4.6 

27.76 

12.50 

9-64 

5O.5 

5-2 

23.8 

4-8 
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TABLE II." ESTIMATED VALUES OP 60 ^/ 6 C L AND 
FESULTANT NEUTRAL POINT 


Airplane 

Wing 

aerodynamic 

center 

(fraction c) 

Estimated values of increment in pitching- 
moment slope, dC^/dCx, 

Airplane neutral 
point 

( fraction - cT) 

V C L/ 

/ 60 A 

W, 

4 c t./L 0 

& 

•'Pi 

w 4 

Q 

Estimated 

*0 

Flight 

x o 

1 

0.245 

0.110 

0.054 


0.006 

0.007 

-O.I77 

O.36 

O.36 

2 

• 2k6 

.084 

.047 


.013 

.013 

-.157 

■ 53 

.32 

3 

. 246 

.094 

.031 


.019 

.012 

-.156 

■34 

•32 

4 

. 23 g 

.113 

.029 


.022 

.011 

-.175 

■35 

•37 

5 

• 2*8 

.103 

.051 


.023 

.019 

-.196 

■ 34 

•33 

6 

.240 

. 077 - 

.027 


.0.19 

.011 

-. 13 k 

■ 32 

• 33 

7 

. 24 ° 

.087 

.021 


.030 

.012 

-.150 

•35 

.32 

8 

,2kl 

.170 

.048 

0.040 

.030 

.025 

-.313 

•41 

.41 

9 

. 24-5 

.111 

.°39 


.011 

.012 

-.173 

•35 

• 33 

10 

.238 

.128 

.048 


.018 

.016 

-.210 

•57 

.36 

11 

.238 

.090 

.039 

— 

.017 

.012 

-.158 

• 33 

• 34 

12 

.238 

.108 

.037 


.013 

.014 

-.172 

•35 

• 34 
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Center of gratify > fraction c 


Figure l.-I/iustration of procedure for determining neufroi point 

for iypfca / a/ rg/ane - 
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NACA CB NO. L5C01 

Fig. £ 
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a) For use tv/fh fuselage section directly ahead of wing 


.8 _ i.2. 


v) ror use with ail fuselage sections ahead of wing except 
section adjacent to wine. - 


Figure 2~ Variation of dS/dac with distance from w/ng 
leading edge . j S- = 4*5 radians. (Data from references) 
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Figure J.- Variation of dC N /dec and dC L /dx. w/ti) aspect ratio. (Dote 
fro/77 reference 4 -.*) 
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Fig. 4 


Figure 
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End- plat e height 
Horizontal- fail span 


4.— Correcfion for dC^jdaij. of horizon 
end p/pfea. (Doha? from reference 4.) 
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Figure S. - 
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Views of airplanes tested. 
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fTsfimated neatrot point, percent c 


NACA CB No. L5C01 


Fig. 6 
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F/ig/jt neutra/ point, percent c 


Figure < 5 . - Comparison of estimated stick- fixed 
neutrai points witt) neutrai points from 
ftiytjf tests. Cj_ - 0.3. (jest-po/nt symbois 
are / dentif/ed in tab/e Z.) 
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